Spinodal-assisted crystallization in polymer melts 
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Recent experiments in some polymer melts quenched below the melting temperature have re- 
ported spinodal kinetics in small-angle X ray scattering before the emergence of crystalline structure. 
To explain these observations we propose that the coupling between density and chain conformation 
induces a liquid-liquid binodal within the equilibrium liquid-crystalline solid coexistence region. 
A simple phenomenological theory is developed to illustrate this idea, and several experimentally 
testable consequences are discussed. Shear is shown to enhance the kinetic role of the hidden binodal. 
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Upon cooling a crystallizable polymer melt a hierar- 
chy of ordered structures emerges [EJ. First there are 
crystalline 'lamellae', comprising regularly packed poly- 
mer chains, each of which is ordered into a specific helical 
conformation. These lamellae interleave with amorphous 
layers to form 'sheaves', which in turn organise to form 
superstructures (e.g. spherulites). These structures may 
be probed by various techniques: e.g. wide-angle X ray 
scattering (WAXS) is sensitive to atomic order within 
lamellae ('Bragg peaks'), while small-angle X ray scatter- 
ing (SAXS) probes lamellae and their stacking. Inspired 
by recent experiments, we propose a model for the earli- 
est stages of ordering in a supercooled polymer melt, and 
discuss several experimentally-testable consequences, in- 
cluding strain and pressure effects. 

In a supercooled simple liquid, the following picture 
M is widely accepted. Nuclei of the lower free energy 
(crystal) phase are constantly formed by thermal fluc- 
tuations. But the cost of creating an interface means 
that only large enough nuclei can grow — the melt is 
metastable. An induction time, rj, elapses before the 
probability of forming such 'critical nuclei' becomes sig- 
nificant. This picture is usually deemed appropriate for 
polymer melts; instead effort is focused on explaining the 
anisotropic shape and growth rate of crystal nuclei M . 

In the 'classical' picture of polymer melt crystalliza- 
tion we expect and observe Bragg peaks in WAXS after 
an induction period r,. SAXS accompanies the WAXS, 
corresponding to interleaved crystal lamellae and amor- 
phous regions ffl. No SAXS is expected during n. How- 
ever, recent experiments have reported SAXS peaks dur- 
ing the induction period and before the emergence of 
Bragg peaks. Initially the SAXS peak intensity grows 
exponentially, and it may be accurately fitted to Cahn- 
Hilliard (CH) theory for spinodal decomposition — the 
spontaneous growth of fluctuations indicative of thermo- 
dynamic instability The peak moves to smaller an- 
gles in time, stopping when Bragg peaks emerge. By 
fitting to CH theory an extrapolated spinodal tempera- 
ture (at which the melt first becomes unstable towards 
local density fluctuations) can be obtained. Spinodal 



kinetics have been reported in different polymer melts: 
poly(ethylene terphthalate) (PET) pi , poly(ether ketone 
ketone) (PEKK) polyethylene (PE) and isotactic 
polypropylene (i-PP) |7|||- Despite these recent results, 
no coherent model exists for these phenomena. 

Such observations can be explained by appealing to 
the concept of a 'metastable phase boundary', a common 
strategy in metallurgy [pUTc]] . Consider an alloy quenched 
into a region of its phase diagram where we expect co- 
existence between two phases, say, a + j3. If the two 
phases are symmetry-unrelated there is always an en- 
ergy barrier for phase separation, and we do not expect 
to see spinodal (unstable) dynamics. Nevertheless, spin- 
odal dynamics and textures have been observed in some 
such cases. Some time ago, Cahn suggested |J that this 
could be due to a metastable phase boundary for coexis- 
tence between two symmetry- related phases (say, a + a') 
buried within the equilibrium a + (3 coexistence region. 
Recently, similar ideas are used to explain slow kinetics 
in the formation of colloidal crystals in colloid-polymer 
mixtures and globular protein crystallization |l2| . 

Similarly, a plausible explanation for the observation 
of spinodal dynamics in polymer melts is the presence of 
a metastable liquid-liquid (LL) phase coexistence curve 
(or 'binodal') buried inside the equilibrium liquid-crystal 
coexistence region, Fig. |l| Quenching sufficiently below 
the equilibrium melting point T m , we may cross the spin- 
odal associated with the buried LL binodal at tempera- 
ture T s < T m . Below, we give a physical mechanism that 
can give rise to such a metastable binodal, calculate the 
phase diagram using a phenomenological free energy, and 
delineate some consequences of our model. 

Our starting point is the unremarkable statement that, 
in order to crystallize, polymer chains must adopt the 
correct conformation. For example, chains in crystalline 
PE have the all-trans (or 'zig-zag') conformation, while 
in the melt the conformation is randomly trans or gauche. 
Generally, the preferred conformation is some form of he- 
lix. Furthermore, the radius of gyration of a (very long) 
chain changes little during crystallization, suggesting |l3| ] 
that neighboring segments adopt the correct conforma- 
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tion and crystallize 'in situ'. It is commonly assumed that 
conformational (intrachain) and crystalline (interchain) 
ordering occur simultaneously. We suggest, in light of 
recent experiments, that these processes can occur se- 
quentially. To motivate this suggestion we examine more 
closely the physics of conformation changes. 



Up, p*) = p [h (A T ) pI + ¥pI + {c P i] 
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FIG. 1. Proposed generic phase diagram for a polymer melt 
calculated as described in the text. T m and T s are the melt- 
ing and spinodal temperatures encountered along the (con- 
stant density) quench path (dotted line). Parameters used 
are RM b = 0.8,fc B T, = 0.29 E ,v = lAE w, A = Q.lauo, 
b = -0A(v a 3 /E o ) 1/2 , c = 0.5 a 2 v /E , and a = 0.8 /w. In- 
set shows the measured induction time as a function of tem- 
perature for isotactic polypropylene [jjj. 

In a melt, it is believed that chain conformation alone 
cannot drive a phase transition. However, conforma- 
tion is coupled to density. Chains with the 'correct' 
(helical) conformation typically pack more densely than 
those with more or less random conformations. More- 
over, the energy barriers between different rotational iso- 
meric states (RIS) are density-dependent 0. We now 
show that conformation- density coupling can induce a LL 
phase transition. A phenomenological free energy which 
incorporates these effects is a function of the following or- 
der parameters: the (average) mass density p; the coeffi- 
cients {p q } in the Fourier expansion of the crystal density 
in terms of the appropriate stars of reciprocal lattice vec- 
tors {q} (essentially the intensities of Bragg peaks) [ff5| ; 
and the occupancies {r/i} of various RIS (and thus chain 
conformation). To illustrate the principles, we assume 
that a single p q = p* and a single r\ suffice, correspond- 
ing to a fictitious polymer with body-centered cubic crys- 
tal structure [|l6| and two RIS. The free energy per unit 
volume has three components: 



/ = fo(p) + f*(p,P*) + fviViPiP*)- 



(1) 



The equation of state pdf/dp—f —p determines the T—p 
isobar at pressure p. The first term, /o, is the free energy 
of a melt with random chain conformations. Equation- 
of-statc fits to polymer liquids suggest the following form 



f (p) = Rk B Tpln[(l/p)-w] 



(2) 



where R and w are widely tabulated |T7| ]. The (bare) 
Landau free energy of crystallization is taken to be [^5|Jl6| 



For simplicity we let a(p,T) = aok B [T — T* (1 + cup)], 
where a and T* account for the enhancement of crystal- 
lization due to increased density, fo + f* has a double well 
structure and gives a (bare) first-order transition between 
amorphous (p* = 0) and crystalline (p* ^ 0) states. 

f v describes how the distribution of chain conforma- 
tions varies smoothly from random (rj = 0) to totally 
ordered (helix, n = 1) as the temperature is lowered to 
zero Jl8[ . In isolation a polymer thermally populates its 
RIS with a Boltzmann distribution P a ~ exp{— j3E a }, 
where E a is the energy of state <f> a relative to the ground 
state <f)Q, and (3 — l/k B T. As the temperature is low- 
ered the mean occupancy r\ of state 4>o, relative to the 
T = 00 disordered state (ry = 0), increases. We describe 
this process, for a two-state model, by 



fv(v,P,P*) 



k B Tp 
2M b 



f cosh 2 ( — ) - nsinh(l3E) 
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where Mb is the mass of a monomer with characteris- 
tic volume vq = wMf,. Minimizing f n over r\ yields the 
correct Boltzmann distribution r](T) = tanh(/3i?/2). We 
choose 



E(p, p*) = Eg + vp + Xpl. 



(5) 



As more bond sequences occupy the ground state, 
monomers can rearrange to pack tighter and reduce the 
excluded volume interaction (hence the perturbation vp). 
A positive v encourages phase separation to take advan- 
tage of this density-conformation coupling. Similarly, 
adjacent ground state sequences enhance crystallization 
(hence the term ApJ). The A term is quadratic in p* by 
symmetry fl6f| . 

To calculate the phase diagram in the temperature- 
densiy (T — p) plane we first minimize / with respect to 
77. Note that minimizing f v over 77 and expanding / in p 
renormalizes a(p, T) and c in Eq. (|J) and, at sufficiently 
small T, destabilizes the homogeneous melt. Physically, 
the system gives up conformational entropy to relieve 
packing frustration, and separates into a dense, more or- 
dered liquid and a less dense and less ordered liquid. The 
renormalization of c lowers the barrier to crystallization 
in the dense, high-77 liquid. Next, / is minimized with re- 
spect to p* to give a final free energy with two branches, 
/(p, p* = 0) (liquid) and /(p, p* 7^ 0) (solid). The com- 
mon tangent constructed between these two branches at 
any temperature gives the densities of coexisting liquid 
and crystal phases at that temperature, Fig. ^ j20). At 
low enough temperatures, the liquid branch gives a spin- 
odal buried entirely within the equilibrium liquid-crystal 
coexistence region given as usual by the points of inflec- 
tion, d 2 f/dp 2 = 0. If a melt is quenched inside this 
spinodal, it will phase separate into two coexisting liq- 
uids, given by the common tangent construction, with a 
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coarsening interconnected domain texture, Fig. ||. The 
two liquids differ in their distributions of conformations, 
with the denser liquid adopting a distribution closer to 
that needed for crystalline packing. 
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FIG. 2. Schematic free energy density plots at a particu- 
lar temperature. At this temperature, a melt with density 
Pl < p < Pc will phase separate into coexisting liquid and 
crystal phases at densities pL and p c . The common tangent 
drawn on the liquid branch alone gives the densities, pLi and 
PL2, on the metastable LL binodal, for which see Fig. |l|. 

To calculate Fig. |l] we fix the dimensionless parameter 
RMf, = 0.8, appropriate for PE jL7|, leaving v/{Eqw) as 
the only adjustable parameter to determine the spinodal 
temperature T s . For v — IAEqw we find an LL critical 
point at k B T c = 0.75Eq and p c w = 0.53. Experiments on 
PE, for which pw ~ 0.685 (w = 0.875 cm 3 g" 1 E% and 
p = 0.783 gem" 3 give k B T s ~ 0.86 E fffusrng 

Eq = 930 cal/mole |^4[). The value for v corresponds to a 
relative change of E by of order +0.58 on going from sin- 
gle chain to melt conditions, Eq. [|. This agrees, in sign 
and magnitude, with the known behaviour of common 
polymers ]23| . The agreement over sign is particularly 
significant. We suggest that density and conformation 
work cooperatively, i.e. v > 0; no LL binodal within 
the required temperature range was obtained for v < 0. 
Even a crude Landau theory, therefore, puts constraints 
on allowable physical mechanisms. Our choice of crys- 
tallization coefficients (a, b, c, A) (caption, Fig. [I]), gives 
a reasonable value for T m , but the crudeness of Landau 
theory for first order transitions renders detailed fitting 
somewhat meaningless. 

One of the coexisting liquids is closer in density and 
conformation to the crystal phase than the original melt, 
and has a lower energy barrier, A(T, p), to crystalliza- 
tion, so inducing 'spinodal-assisted nucleation'. We ex- 
pect A(T,p) to decrease with increasing quench depth 
below T s (and hence widening LL coexistence gap). In- 
deed, in our simplified model, we find A(T, p) = at 
temperature T p <T S . The induction time, Tj, is then 



const, x e 



A(T,p)/k B T 
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existing LL densities. We expect r s to be weakly (at most 
power law) dependent on temperature. The exponential 
term accounts for the barrier to nucleate a crystal from 
the dense liquid. The strong temperature dependence of 
Ti should change over to a much weaker dependence at 
some T p < T s , where A(Tn,p) < k B T p . This has been 
found in i-PP (inset, Fig. g) §. 

The characteristic length scale associated with the de- 
veloping spinodal texture gives rise to a SAXS peak, 
which evolves initially according to CH theory . The 
coarsening of this texture is observed to be arrested at 
the end of the induction period (typical scale £ m ), when 
Bragg peaks appear in WAXS . It is at present unclear 
how the spinodal texture at the end of the induction pe- 
riod evolves into spherulites. However, the final spinodal 
texture length scale £ m evidently controls the thickness 
of the first crystal lamellae. Moreover, large stress will 
develop once one of the two liquids in a bicontinuous 
texture, Fig. ||, crystallizes. We expect such a texture to 
fragment into individual crystalline lamellae. 
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where r s is the time to reach an intermediate spinodal 
texture, Fig. ||, in which the regions have (almost) the co- 



Amorphous 

FIG. 3. Schematic representation of the late-stage spinodal 
texture for coexisting liquid phases with different conforma- 
tions, showing a single chain; thin line = disordered confor- 
mation, thick line = correct (helical) conformation for crys- 
tallization. Each chain is a 'conformational copolymer'. 

Our arguments so far have been based on 
conformation-density coupling. Once a polymer seg- 
ment has adopted the correct (helical) conformation its 
persistence length should increase, which couples to the 
orientational order of chains. Indeed, depolarized light 
scattering by Imai and coworkers has suggested the ex- 
istence of orientational fluctuations during the spinodal 
phase of a crystallizing PET melt [p2| . Provided that 
orientational ordering is not strong enough to induce a 
separate transition, then the inclusion of a nematic or- 
der parameter in Eq. |l| only renormalises the coefficients 
in fy-dependent terms. In some cases, the increasing 
chain stiffness accompanying conformational order may 
be sufficient to drive an isotropic — > nematic transition, 
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resulting in a three-step process: melt — > (isotropic) liq- 
uid(l) + liquid(2), followed by liquid(2) — > nematic — > 
crystal. This possibility should be investigated. It should 
also be possible to model specific polymers, e.g. using 
density functional theor y p a ] , augmented to include the 
effects leading to Eq. || [|l4|]. 

Until recently, spinodal scattering was mainly observed 
in polymer melts crystallizing under shear p\pjj^^. This 
may be understood in a natural way within the present 
framework. Shear (and extensional) flow couples prin- 
cipally to the orientation of polymer segments, hence 
straightening chains and enhancing r\, thereby biasing the 
tendency towards LL separation. A simple way to incor- 
porate this is to renormalise the activation energy E as 
E — vq<j where a is the stress. It is highly suggestive 
that, for appropriate values of stress under strong flow 
(the plateau modulus Go) and volume (vq above), the 
LL binodal of Fig. I is shifted upward significantly (by 
ST S ~ 0.01Eo/k B ). Flow will shift the liquid-solid coexis- 
tence curve much less because the regions with crystalline 
order will resist deformation. 

Our simple theory suggests several interesting experi- 
ments. First, and most directly, conformational fluctua- 
tions just above T s could be detected and studied, e.g., by 
Raman spectroscopy jjj , perhaps simultaneously with de- 
polarised light scattering (to monitor orientational fluc- 
tuations). Secondly, on approaching a spinodal line, var- 
ious properties (e.g. correlation length) should exhibit 
power-law divergences. The observation of such diver- 
gences will lend much weight to the correctness of our 
model. Thirdly, the point at which the LL spinodal is 
encountered in a quench can be modified by pressure, 
Fig. [T| In particular, it may be possible to access the 
LL critical point, T c : recent simulations suggest a mas- 
sive enhancement of the nucleation rate in the vicinity 
of T c Q. More generally, the coupling of density to 
(molecular) structural order parameters is an emerging 
generic theme in the study of supercooled liquids (wa- 
ter — > amorphous ice, |p6[ ; polymer melts near the glass 
transition [27]]). Finally, processes such as surface nucle- 
ation could give rise to a SAXS peak, but are unlikely 
to follow Cahn-Hilliard kinetics. Also, effects other than 
conformation and orientation (e.g. polydispersity) may 
induce LL phase separation. Experiments on monodis- 
perse alkanes are under way to address this possibility. 
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